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Recently, cold atomic Fermi gases with the large magnetic dipolar interaction have been laser cooled down to 
quantum degeneracy. Different from electric-dipoles which are classic vectors, atomic magnetic dipoles are 
quantum-mechanical matrix operators proportional to the hyperfine-spin of atoms, thus provide rich 
opportunities to investigate exotic many-body physics. Furthermore, unlike anisotropic electric dipolar 
gases, unpolarized magnetic dipolar systems are isotropic under simultaneous spin-orbit rotation. These 
features give rise to a robust mechanism for a novel pairing symmetry: orbital p-wave (L = 1) spin triplet (S 
= 1) pairing with total angular momentum of the Cooper pair / = 1. This pairing is markedly different from 
both the 3 He-B phase in which / = 0 and the 3 He- A phase in which /is not conserved. It is also different from 
the p-wave pairing in the single-component electric dipolar systems in which the spin degree of freedom is 
frozen. 

Ultracold atomic and molecular systems with electric and magnetic dipolar interactions have become the 
research focus in cold atom physics 1 " 7 . When dipole moments are aligned by external fields, dipolar 
interactions exhibit the d r 2_ 3z 2-type anisotropy. The anisotropic Bose-Einstein condensations of dipolar 
bosons (e.g. 52 Cr) have been observed 8 " 11 . For the fermionic electric dipolar systems, 40 K- 87 Rb has been cooled 
down to nearly quantum-degeneracy 3 . Effects of the anisotropic electric dipolar interaction on the fermion many- 
body physics have been extensively investigated. In the Fermi liquid theory, both the single particle properties and 
collective excitations exhibit the d r i_ 3z 2 anisotropy 12 " 17 . In the single -component Fermi systems, the leading order 
Cooper pairing instability lies in the p-wave channel, which is the simplest one allowed by Pauli's exclusion 
principle. The anisotropy of the electric dipolar interaction selects the instability in the p z - channel, which is 
slightly hybridized with other odd partial wave channels 20 " 27 . For two-component cases, the dipolar interaction 
leads to anisotropic spin -triplet pairing, and its orbital partial wave is again in the p z - channel 28 " 31 . The triplet 
pairing competes with the singlet pairing in the hybridized s + d r 2_ 3z 2- channel. The mixing between the singlet 

K 

and triplet pairings has a relative phase + -, which leads to a novel time- reversal symmetry breaking Cooper 
pairing state 29 . 

An important recent experimental progress is the laser cooling and trapping of magnetic dipolar fermions of 
161 Dy and 163 Dy with large atomic magnetic moments (lO^) 1 ' 2 . There are important differences between magnetic 
and electric dipolar interactions. Electric dipole moments are essentially non-quantized classic vectors from the 
mixing between different rotational eigenstates with opposite parities, which are induced by external electric 
fields 3,4 , thus electric dipoles are frozen. In the absence of external fields, even though at each instant of time there 
is a dipole moment of the heteronuclear molecule, it is averaged to zero at a long time scale. In contrast, magnetic 
dipole moments of atoms are intrinsic, proportional to their hyper-fine spins with a Lande factor. Unpolarized 
magnetic dipolar Fermi systems are available, in which dipoles are defrozen as non- commutative quantum 
mechanical operators, thus lead to richer quantum spin physics of dipolar interactions. Furthermore, the mag- 
netic dipolar interaction is actually isotropic in the unpolarized systems. It is invariant under simultaneous spin- 
orbit rotations but not separate spin or orbit rotations. This spin-orbit coupling is different from usual single 
particle one, but an interaction effect. It plays an important role in the Fermi liquid properties such as the 
unconventional magnetic states and ferro-nematic states predicted by Fregoso et a/ 1819 . 

It is natural to expect that magnetic dipolar interaction brings novel pairing symmetries not studied in 
condensed matter systems before. The systems of 161 Dy and 163 Dy are with a very large hyperfme spin of 

F = — , thus their Cooper pairing problem is expected to be very challenging. As a first step, we study the simplest 
2 1 

case of spin--, and find that the magnetic dipolar interaction provides a novel and robust mechanism to the p- 
wave (1=1) spin triplet (S = 1) Cooper pairing to the first order of interaction strength, which comes from the 
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attractive part of the magnetic dipolar interaction. In comparison, 
the p- wave triplet pairing in usual condensed matter systems, such as 
3 He 32 " 34 , is due to the spin-fluctuation mechanism, which is at the 
second order of interaction strength (see Refs. 35,36 for reviews). This 
mechanism is based on strong ferromagnetic tendency from the 
repulsive part of the 3 He- 3 He interactions. Furthermore, the p-wave 
triplet Cooper pairing symmetry patterns in magnetic dipolar sys- 
tems are novel, which do not appear in 3 He. The orbital and spin 
angular momenta of the Cooper pair are entangled into the total 
angular momentum / = 1, which is denoted as the /-triplet channel 
below. In contrast, in the 3 He-B phase 33 , L and S are combined into / 
= 0; and in the 3 He-A phase, L and S are decoupled and /is not well- 
defined 32,34 . There are two competing pairing possibilities in this /- 
triplet channel with different values of J z : the helical polar state (J z = 
0) preserving time reversal (TR) symmetry, and the axial state (J z = 
± 1) breaking TR symmetry. The helical polar state has point nodes 
and gapless Dirac spectra, which is a time- reversal invariant general- 
ization of the 3 He-A phase with entangled spin and orbital degrees of 
freedom. In addition to usual phonon modes, its Goldstone modes 
contain the total angular momentum wave as entangled spin-orbital 
modes. 

Results j 

We begin with the magnetic dipolar interaction between spin— fer- 
mions ~ 

^{Soo''%-3(S aa r-f)(S^-f)} 5 (1) 

where r is the relative displacement vector between two fermions; fi is 
the magnitude of the magnetic moment. Such an interaction is 
invariant under the combined SU(2) spin rotation and SO(3) space 
rotation. In other words, orbital angular momentum L and spin S are 
not separately conserved, but the total angular momentum J = L + S 
remains conserved. Its Fourier transformation reads 19 

V*ptfAq) = y ^{^(Saa'-q) {Stf-q) (2) 
The Hamiltonian in the second quantization form is written as 

1 (3) 

+ 2V E VxM '«' (*~ (* ; ? ) Pp «' ? ) ' 

k,k',q 

where e(]cj =hk 2 /(2m); fi c is the chemical potential; 

Pp> a , (k\ c^J =cp ^ — k + cfjca' (k + q^j is the pairing operator; the 

Greek indices a, a' and refer to j and [; V is the volume of 
the system. We define a dimensionless parameter characterizing the 
interaction strength as the ratio between the characteristic inter - 

. 2 ii 2 mkf 

action energy and the Fermi energy: A = E int E F = —A 

3 n 2 h 

We next study the symmetry of the Cooper pairing in the presence 
of Fermi surface, i.e., in the weak coupling theory. An important 
feature of the magnetic dipolar interaction in Eq. (1) is that it 
vanishes in the total spin singlet channel. Thus, we only need to study 
the triplet pairing in odd orbital partial wave channels. Considering 
uniform pairing states at the mean-field level, we set q = 0 in Eq. (3), 

and define triplet pairing operators P s (itj , which are eigen-operators 

of Siz + S 2z with eigenvalues s z = 0, ± 1, respectively. More explicitly, 

they are P 0 (k) = ±= [p n (k) +P it (k)] , P 1 (k)=P^(k), 

P_i (itj =P[[ (icj . The pairing interaction of Eq. (3) reduces to 



k,k',s z s' 



(4) 



where 



^(^'HE^ 



(5) 



l5 2 



1 1 



1 



a - p ) is the Clebsch-Gordan coefficient for two spin-- states 



to form the spin triplet; and V SzS > z (k\ k'^J is an odd function of both k 
and k! . 

The decoupled mean-field Hamiltonian reads 

where we only sum over half of the momentum space; 
^{k^=e(^j— fi ch and \i ch is the chemical potential; 

^P^ifc) = (c T (ifc) ^ifc) ( — ifc),cj( — A aP is defined as 



tion as 



^ oc^P ) * A 5z . A 5z satisfies the mean-field gap func- 



(7) 



(7) 



k',s' z 

where K (k'^j = tanh ^ E { (k'^j j \lE { The integral in Eq. 

is already normalized following the standard procedure 20 . For sim- 
plicity, we use the Born approximation in Eq. (7) by employing the 
bare interaction potential rather than the fully renormalized T- 
matrix, which applies in the dilute limit of weak interactions. The 
pairing symmetry, on which we are interested below, does not 
depend on the details that how the integral of Eq. (7) is regularized 
in momentum space. The Bogoliubov quasiparticle spectra become 

E\,2 {k^j = + 1\ 2 {k^j , where k\ 2 (k^j are the eigenvalues of the 

positive-definite Hermitian matrix A^ (k^j A (k^j . The free energy can 
be calculated as 



F = — — In 2 cosh 

/Lf* I 2 



k,i= 1,2 



(8) 



2V^ 



where vj^ 1 (k; k'^j is the inverse of the interaction matrix defined as 

k',s' z 

We next linearize Eq. (7) around T c and perform the partial wave 
analysis to determine the dominant pairing channel. Since the total 
angular momentum is conserved, we can use / to classify the 
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eigen-gap functions denoted as <j> a s f z (k^J . The index a is used to 

distinguish different channels sharing the same value of /. (j) a s f z (k^j 
satisfies 

N ° J H ( l > *T (?) = ^ (*) • ( 10 ) 

where N 0 = — is the density of state at the Fermi surface; wf are 

n 2 h J 

dimensionless eigenvalues; k, k! are at the Fermi surface. Then Eq. (7) 
is linearized into a set of decoupled equations 

(11) 



f Jh { 1 + w a j [ln(2^Q) / (7dk B T)] } = 0, 



where Q is an energy scale at the order of the Fermi energy playing 
the role of energy cut-off from the Fermi surface. 

The decomposition of Vj^ (k\ k!^j into spherical harmonics can be 
formulated as 

VL mS ,,L>rr,wY* Lm (n k )Y Um ,(n i ), (12) 

Lm,L'm' 



N 0 T 
— V T , 
All 



where L = V orL = L'±2, and L, V are odd numbers. The expres- 
sions of the dimensionless matrix elements V LmSz . L > m > s > z are lengthy 
and will be presented elsewhere. By diagonalizing this matrix, we find 
that the most negative eigenvalues is w J=l = —3nX/4 lying in the 
channel with / = L = 1. All other negative eigenvalues are signifi- 
cantly smaller. Therefore, dominate pairing symmetry is identified as 
the /-triplet channel with L = S = 1 in the weak coupling theory. 
Following the standard method in Ref. 20 , the transition temperature 
2e y Cl l — 

T c is expressed as T c ~ e I wI ~ 1 1 . For a rough estimation of the order 

K 

of magnitude of T a we set the prefactor in the expression of T c as Ef. 

In order to understand why the /-triplet channel is selected by the 
magnetic dipolar interaction, we present a heuristic picture based on 
a two-body pairing problem in real space. Dipolar interaction has a 
characteristic length scale a dp = mpflti 2 at which the kinetic energy 
scale equals the interaction energy scale. We are not interested in 
solving the radial equation but focus on the symmetry properties of 
the angular solution, thus, the distance between two spins is taken 
fixed at a^p. We consider the lowest partial-wave, p-wave, channel 
with L= 1. The 3X3 = 9 states (L=S= 1) are classified into three 
sectors of /= 0, 1 and 2. In each channel of /, the interaction energies 
are diagonalized as 

1 1 



Eo=Edp, Ei = — -Edp, E 2 = —Edp, 



(13) 



respectively, where E^p = j a\ p . Only the total angular momentum 

triplet sector with /= 1 supports bound states, thus is the dominant 
pairing channel and is consistent with the pairing symmetry in the 
weak- coupling theory. 

This two -body picture applies in the strong coupling limit. 
Although a complete study of the strong coupling problem is beyond 
the scope of this paper, this result provides an intuitive picture to 
understand pairing symmetry in the /-triplet sector from spin 

configurations. We define that ^ and pp (&J are eigenstates with 

eigenvalues zero for operators ^-(5i+S 2 ) and e fl -L(fi = x,y,z) 1 
which are the total spin and orbital angular momenta projected 

along the e^-direction. The /-triplet sector states are <j> fl (Q) = 

-^ e pva v P*(ty with K satisfying (^-7)^ = 0. For example, 



=^™ fl {k>ik>2+|^) 1 K) 2 }> 



(14) 




Figure 1 | The spin configurations of the two-body states with a) / = 1 
and j z = 0 and b) / = j z = 0. The interactions are attractive in a) but 
repulsive in b). 



where e p =£cos</>+j>sin^ and \a 6p ) and \Pe p / are eigenstates of 

e p -a with eigenvalues of ±1. As depicted in Fig. 1 A, along the 
equator where <j) z has the largest weight, two spins are parallel and 
along r, thus the interaction is dominated by attraction. On the other 
hand, the eigenstate of / = 0 reads 

*o(«)=x^(n) = ^{l^>ilfo> 2 + l&)il^)2}, ( 15 ) 

where |a Q ) and |jS Q ) are eigenstates of Q a with eigenvalues ±1. As 
shown in Fig. 1 B, along any direction of Q, two spins are anti-parallel 
and longitudinal, thus the interaction is repulsive. 

Let us come back to momentum space and study the competition 
between three paring branches in the /-triplet channel under the 
Ginzburg- Landau (GL) framework. We define 

M fc >;7f[- Al (*) +A ->(*)]' 

(16) 



A ' (*) = ^ [ Al (*) + A " 1 (*) ] ' Az (*) = Ao (*) • 

The bulk pairing order parameters are defined as = ^ 

kp (k)j , where no summation over \x is assumed. We define pairing 
parameters and their real and imaginary parts as the following 3- 
vectors A = (A X ,A 7 ,A Z ) . The GL free energy is constructed to main- 
tain the (7(1) and SO(3) rotational symmetry as 

F = aA*-A + y 1 |A*-A| 2 + y 2 |A* x a| 2 , (17) 



where 



a = iV 0 In 



(18) 



The sign of y 2 determines two different pairing structures: ReA||lmA 

at y 2 > 0, and ReA_LImA at y 2 < 0, respectively. Using the analogy of 
the spinor condensation of spin-1 bosons, the former is the polar 
pairing state and the latter is the axial pairing state 37 " 40 . 

For the polar pairing state, the order parameter configuration can 
be conveniently denoted as A = \ A \z up to a U( 1 ) phase and SO(3) - 
rotation. This pairing carries the quantum number J z = 0. The 

pairing matrix A p J p = X - |A| \k y O\-k x o 2 )io 2 \^ reads 



:i|A| 

2 1 1 



[k y + ik^j 



ky ik x 



(19) 
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It equivalents to a superposition ofp x + ip y orbital configurations for 
spin- 1 1 (||) pairs, respectively. Thus, this pairing state is helical. It is 
a unitary pairing state because A is proportional to a 2X2 identity 
matrix. The Bogoliubov quasiparticle spectra are degenerate for two 

different spin configurations as = ^/^+ | a" (it) | with the 

anisotropic gap function |A p/ ^| 

Fig. 2. They exhibit Dirac cones at north and south poles with oppos- 
ite chiralities for two spin configurations. 

Similarly, the order parameter configuration in the axial pairing 

z^\A\(e x -\-iey) up to the symmetry 

l.The 



^ |A| z sin 2 0fc depicted in 



state can be chosen as A = 



V2 



transformation. This state carries the quantum number of J 2 
pairing matrix A^ = — - | A 1 { \k z {o x + io 2 ) + a z (k x + ifcyj io 2 } 



takes the form 



-ih 



2/i 



(20) 



This 
|A| 



a non-unitary pairing state since 



A'A = 



tions as E 

|2 



. The Bogoliubov quasiparticle spectra 
have two non- degenerate branches with anisotropic dispersion rela- 
™ 2 (fj = y^+|A fl ^^| 2 . The angular gap distribution 

I A ± (^) | = ~ I A ! 2 ( : ± cos ^) 2 is depicted in Fi g- 2 - Eacn of branch 
1 and 2 exhibits one node at north pole and south pole, respectively. 
Around the nodal region, the dispersion simplifies into 

Ei, 2 (jc) = y ^ vj(k z + kf) 2 + ^ \A\ 2 (k\\/kf) 4 , which is quadratic in 



the transverse momentum k\\ 



k%-\-ky. 



At the mean-field level, the helical polar pairing state is more stable 
than the axial state. Actually, this conclusion is not so obvious as in 
the case of 3 He-B phase, where the isotropic gap function is the most 
stable among all the possible gap functions 33 . Here, the gap functions 
are anisotropic in both the polar and helical pairing phases. We need 
to compare them by calculating their free energies in Eq. (8). The 
second term contributes the same to both pairing phases. Thus, the 
first term determines the difference in free energies. Let us define the 
ratio between angular integrals of the free energy kernels in Eq. (8) of 
the two phases as 



y(h, k 2 )-- 



JdQ* 2 In 


2 cosh - \ 
2 V 


e k +\A pi (k) 


2 


JdQ, £ ± ln 


2 cosh^i 




A^(fc) 2 



where X\ = 



1 



1 



(21) 



r.y(ki, X 2 ) is numerically plotted in Fig. 3. 



For arbitrary values of £ k , and |A|, y is always larger than 1. 
Therefore, the polar state is favored more than the axial state. This 
can be understood from the convexity of the nonlinear term in 

Eq. (8), which favors isotropic angular distributions of | 42 • 

Although neither gap function of these two states is absolutely iso- 
tropic as in the 3 He-B phase, the polar gap function is more isotropic 
from Fig. 2 and thus is favored. However, we need to bear in mind 
that we cannot rule out the possibility that certain strong coupling 
effects can stabilize the axial state. In fact, the 3 He-A phase can be 
stabilized under the spin feedback mechanism 35 , which is a higher 
order effect in terms of interaction strength. 

Next we discuss the classification of Goldstone modes and vortices 
in these two states. In the helical polar state, the remaining symmet- 
ries are SOj(2) X Z 2 as well as parity and time-reversal (TR), where Z 2 
means the combined operation of rotation n around any axis in the 
xy-plane and a flip of the pairing phase by n. The Goldstone manifold 
is 

[SO/(3) x U c (l)]/[SOj(2)®Z 2 ] = [Sj x U c (l)]/Z 2 . (22) 

The Goldstone modes include the phase phonon mode and two 
branches of spin-orbital modes. Vortices in this phase can be clas- 
sified into the usual integer vortices in the phase sector and half- 
quantum vortices combined with 7i-disclination of the orientation of 
A. In the axial state, the rotation around z-axis generates a shift of the 
pairing phase, which can be canceled by a U c (l) transformation, thus, 
the remaining symmetry is SO; z _^(2). The Goldstone manifold is 
S 2 X U c (l). Only integer vortices exist. 

Discussions 

In summary, we have found that the magnetic dipolar interaction 
provides a robust mechanism at first order in the interaction strength 
for a novel p- wave (1=1) spin triplet (S = 1) Cooper pairing state, in 
which the total angular momentum of the Cooper pair is / = 1 . This is 
a novel pairing pattern which does not appear in 3 He, and, to our 
knowledge, neither in any other condensed matter systems. These 
pairing states include the TR invariant helical polar pairing state and 




cos8k 

Figure 2 | The angular distribution of the gap function | A {k^j | v.s. cos0^ 
in the helical polar pairing state (the red line) and the axial pairing state 
(the black line). 




Figure 3 | The ratio of the angular integrals of the free energy kernels 

^ G?j^| ' A^f) ' W ^ C k * s a ^ wa y s l ar g er than 1. This means that the polar 
pairing is favored at the mean-field level. 
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the TR breaking axial pairing state, both of which are distinct from 
the familiar 3 He-A and B phases. 

Many interesting questions are open for further exploration, 
including the topological properties of these pairing states, vortices, 
spin textures, and spectra of collective excitations. The above theory 

only applies for spin-- systems, in which the magnetic dipolar inter- 
action is too small. For the pairing symmetry in a magnetic dipolar 
system with a large spin S, our preliminary results show that the basic 
features of the /-triplet pairing remains. The spins of two fermions 
are parallel forming S tot = 2S with orbital partial- wave L = 1, and the 
total / = 2S. In the current experiments in Ref. 41 , the highest attain - 

21 

able density reaches 4X 10 13 cm 1 for 161 Dy atoms with S = — . The 

corresponding dipolar energy is E int ~ 2 nK and the Fermi energy for 
unpolarized gases 13.6 nK, and thus X = E int /Ef~ 0.15. If we use 
the same formula of w ]=l above for an estimation of the most negative 
eigenvalue, we arrive at T c /7y~ 0.06, which means that T c ~ 0.8 nK. 
Although it is still slightly below the lower limit of the accessible 
temperature in current experiments, we expect that further increase 
of fermions density, say, in optical lattices will greatly increase T c . 

Method 

We have used the methods of the symmetry analysis, strong coupling analysis, mean- 
field theory, partial- wave analysis, and the Ginzburg- Landau free energy, which have 
been explained in Sec. I. 
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